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ABSTRACT 

Galaxy cluster surveys will be a powerful probe of dark energy. At the same time, cluster abundance 
is sensitive to any non-Gaussianity of the primordial density field. It is therefore possible that non- 
Gaussian initial conditions might be misinterpreted as a sign of dark energy or at least degrade the 
expected constraints on dark energy parameters. To address this issue, we perform a likelihood analysis 
of an ideal cluster survey similar in size and depth to the upcoming South Pole Telescope/Dark Energy 
Survey (SPT-DES). We analyze a model in which the strength of the non-Gaussianity is parameterized 
by the constant /atl; this model has been used extensively to derive Cosmic Microwave Background 
(CMB) anisotropy constraints on non-Gaussianity, allowing us to make contact with those works. We 
find that the constraining power of the cluster survey on dark energy observables is not significantly 
diminished by non-Gaussianity provided that cluster redshift information is included in the analysis. 
We also find that even an ideal cluster survey is unlikely to improve significantly current and future 
CMB constraints on non-Gaussianity. However, when all systematics are under control, it could 
constitute a valuable cross check to CMB observations. 
Subject headings: cosmology: theory - galaxies: clusters - dark energy 



1. INTRODUCTION 

Of the many fascinating discoveries in cosmology over 
the last decade, perhaps none have aroused more inter- 
est than the discovery of the accelerated expansion o f 
the Universe (Ries s et al.lll998t iPerlmutter et al.lll999D . 
Probing the nature of the dark energy thought to be 
driving this acceleration has become a top priority for 
the community, and among the promising tools under 
consideration are surveys of galaxy clusters. Since the 
number of clusters as a function of redshift and mass 
depends on both the growth of structure and on the 
volume of space, the cluster abundance is sensitive to 
the matter density, the density fluctuation amplitude, 
and the expansion history of the Universe. For this rea- 
son, upcoming cluster surveys will be powerful probes of 
cosmology (e.g. Haiman et al. 2000; Holder et al. 2001; 
Eattve & Wellcr 2003; Molnar ct al. 2004; Wan get all 
12004^ .Rapetti et al. 2005: Ma,ria,n Bernstein. 200^ 

Although constraining dark energy is a leading 
motivator for much of the interest in cluster sur- 
veys, it is worth noting that the cluster abundance 
is potentially sensitive to various cosmological pa- 
rameters beyond those associated with dark energy. 
For example, it has been recognized for some time 
that slight deviations from Gaussianity in the pri- 
mordial matter distribution would cause a significant 
chang e in the high mass tail o f the halo distribu- 
tion llLucchin fc Matarrese] Il988t iColafrancesco et al.l 



IChhi et alJ JT99S 
m & Bakerl 120001 



19891 

Robinsonir^ak er[ 

Robinson ct all l2000t iMatarrese et all l2000|) . In 



IRobinson et alJ 
iKovama et al.l 



1996 
199£ 



Electronic address: |emllla no@f nal . gov| 

^ Particle Astrophysics Center, Fermi National Accelerator Lab- 
oratory, Batavia, IL 60510-0500 

^ William I. Fine Theoretical Physics Institute, University of 
Minnesota, Minneapolis, MN 55455 

^ Kavli Institute for Cosmological Physics and Department of 
Astronomy &; Astrophysics, The University of Chicago, Chicago, 
IL 60637 



this paper, we use a maximum likelihood analysis to 
investigate the extent to which dark energy constraints 
from cluster surveys are degraded by including the pos- 
sibility of non-Gaussian initial conditions, in particular 
when considered within the limits allowed by present 
and future CMB observations. 

The specific form of non-Gaussian initial conditions we 
consider here is of the local type, described in position 
space by a primordial curvature perturbatioii of the form 
(|Verde et al.ll2nont iKomatsu fc SDerEfe]ll200lD 



$(x) 



+ /^i[</.2(x)-(02(x))] 



(1) 



where (/>(x) is a Gaussian random field and the de- 
gree of non-Gaussianity is parameterized in terms of 
the constant //vl- For this model, tight constraints of 
the order of IS. fur, ~ 40 are provided by CMB ob- 
servations llKomatsu et al.ll200a iCreminelh etlill 120061: 



ISnergel et all 120061: iChen fc Szanudil l200fiH . while con- 
straints that are somewhat weaker but that are closer 
in physical scale to that of clusters are ex pected 
from higher-order galaxy correlations ( Scoccimarr o et alJ 
From a theoretical point of view, the 
non-Gaussian model of Eq. |^ is motivated in 
part by studies of the generation of density per- 
turbations in inflationary scenarios; while single- 
fleld inflation models typically pr edict an unob- 
servab ly small va lue for f^L (e.g. lAcauaviva et alJ 
l2003t IMaldacenal l200l . multi-field inflation mod- 
els can lead to much higher values (e.g. iLvth et all 
•2003^ 'Dvali et al.' "200?; "Zaldarriaea "2001: 'CreminellJ 
2003; Arkani-Hamed et al. 200 4; Ahshahiha et al. 200j 
iKolb et al. 2006; S asaki et ai.ll2006f) . For a review, see 
iBartolo et al.. (,,20041) . 

While we believe it is worthwhile to keep an open mind 
to other forms of non-Gaussianity which may not be 
properly described by the simple expression in Eq. 
and which might make the extrapolation of current CMB 
constraints to cluster s cales less stra-igh tforward than we 
assume here (see, e.g.. lMathis et aLll2004f) . we note that 
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the physical scale probed by clusters differs from that of 
the Planck survey by roughly a factor of two, so that the 
two probes are likely to be affected more or less equally 
by deviations from Eq. (|T|l. 

As we discuss below in greater detail, the parameters 
in our likelihood analysis include //vl, the matter den- 
sity flm and the matter fluctuation amplitude erg, while 
we consider both a constant and time-varying dark en- 
ergy equations of st ate described in terms of one (w) 
and t wo [wq and w„ . iChevallier fc Po larski"2 001tlLindeil 
|2003|) parameters respectively. For definiteness, we as- 
sume a fiducial ideal survey similar in size and depth to 
that of the upcoming South Pole Tel escope /Dark Energ y 
Survey (SPT-DES, Ruh l et alJ2dollAbbott et al.l20nfSfl . 
We assume a ACDM fiducial cosmology, for two values of 
(Tg, since cluster number counts are extremely sensitive 
to this parameter. 

This paper is organized as follows. In section |21 we in- 
troduce our model for the non-Gaussian mass function 
and describe our analysis of the dependence of the ex- 
pected errors on cosmological parameters on the non- 
Gaussian component. In section|2|we present our results 
and we conclude in sectional 

2. THE MODEL 

In this section we present the methods applied in the 
present work. We begin with a brief review of previ- 
ous works dealing with non-Gaussian initial conditions 
in galaxy cluster observations, and then we describe in 
detail our treatment of the non-Gaussian mass function. 
We conclude this section with a discussion of the likeli- 
hood analysis whose results will be given in sectional 

2.1. Historical overview 

Expressions for the cluster mass function in the pres- 
ence of non-Gaussian initial conditions have been de- 
rived as extensions to th e Press-Sche chter ansantz (PS, 
iPress fc Sc hcchtcr 197 41) first b y Lucchin fc Matarresd 
1)1988(1 and .Colafrancesco et alJ 11989) while a s impler 
app roach has been adopte d later bv Ichiu et aD l|1998() 
and lRobinson et all l(1998ri . 

The original PS formula describes the comoving num- 
ber density n{M)dM of clusters with mass in the interval 
(M, M + dM) as 



nps{M)dM 



2p d 
MdM 



PGiy)dy 



dM, (2) 



where we suppress, for clarity, the redshift dependence, 
p is the comoving mass density, (Tj\/ is the r.m.s. of 
mass fluctuations in spheres of radius R = {SM/AnpY^^ , 
6c = 1.686 is the critical linear overdensity in the spher- 
ical collapse model and Pq is the Gaussian probabil- 
ity distribution function (PDF), Poiy) = e-y^/'^/V^^- 
Since the function Pciv) does not depend explicitly on 
the mass M, and therefore on the scale R, Eq. (jSJ reduces 
to 

The PS formalism assumes that the scale depen- 
dence of the PDF of the density field is completely de- 
scribed by the scale dependence of the variance aj^j. 



iLucchin fc MatarreseHl988l).lColafrancesco et alJ l|1989n 
and. later lMatarrese et alJ 1(200(1) considered a derivation 
of the non-Gaussian mass function, based on Eq. (|2Jl, 
that takes into account the scale dependence of higher 
order cumulants, thereby allowing for a generic depen- 
dence of the PDF on t he sm oothing scale R. Specif- 
ically, iMatarrese et all l(2000() (hereafter MVJ) derived 
the mass function corresponding to the model described 
by Eq. . The non-Gaussianity of the mass function is 
described, in first approximation, in terms of the skew- 
ness 5*3, ij of the smoothed density field 5r, 
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(4) 



and it is obtained from the cumulant generator of the 
distribution as 



nMVj{M)dM ~ X 



1 51 dS- 



3,R 



din (Tm 



2(5, d\nM 



d\nM 



2tt 



-dM, 



(5) 



where (5* = Sc^Jl - S3^rSc/3. 

It's worth noticing here that although Eq. Q should 
be seen as a truncated expansion in powers of </>, the 
mass function provided by Eq. |SJ) is not linear in the 
non-Gaussian parameter //v^ (since S^.r ~ //vl); rather 
it describes the non-Gaussian PDF by its proper depen- 
dence on the skewness while neglecting all higher order 
cumulants. 

The simpler exte nsion to non-Gaussian initial condi- 
tions introduced bv lGhiu et a l. ( 1998) consists instead of 
replacing the Gaussian function Pciv) in Eq- ® by the 
appropriate, non-Gaussian PDF P^cd/), assumed to be 
scale- independent. The resulting mass function, which 
we will denote here as "extendcd-PS" or EPS, therefore 
reads 



nEPs{M)dM 



M2 CTM 



^ ji . r Png{5c (JM)dM. 
d m M 



(6) 



This approach, has been tested in N-body simulations 
by Robinson fc Baker ( 2000) for several non-Gaussian 
models; they find that Eq. (jHJ agrees with measurement 
of the cumulative mass function n{> M) in the simula- 
tions to within 25%. While this error is slightly larger 
than the differences between the PS formula, Eq. (|2J|, 
and simulation results for Gaussian initial conditions, it 
is much smaller than the model-to-model differences be- 
tween the cumulative mass functions. As a measure of 
the non-Gaus sianity of the ta i l of th e distribution func- 
tion P/vg(2/), ' Robinson et al.l l(1998() introduced the pa- 
rameter G (there called T) defined as 



G 



PNGiy)dy 
irPGiy)dy 



(7) 



with G = 1 corresponding to t he Gaussian case . 

Following this approach 'Robinson et al.l l(1998() . 
Kovama ct al. (1993) and ^illick (200^ placed con- 
straints on primordial non-G aussianity from X-ray 
cluster survey observations (iHenrv fc ArnaudI I199H 
lEbehng et al.il996l:lHenrvlll997() and lAmara fc RefregieJ 
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(|2004(1 relate primordial non-Gaussianity with the nor- 
malization of the dark matter power spectrum. In par- 
ticular, assuming that the non-Gaussian primordial field 
can be Kenerical l y desc ribed by a log-normal distribution, 
iR.obinson et alJ 1)200(11) found, for a ACDM cosmology, 
the constraint G < 6 at 2cr level. An analysis of the con- 
straining power of future Sunyaev-Zel'dovich (SZ) cluster 
surveys on cosmological parameters which includes the 
possibility of primo rdial non-Gaussianity is provided by 
iBenson et al.l l)2002f) . Specificall y, this work a ssumes the 
log-normal PDF studied by Robinson et alJ (f2000j and 
performs a Fisher-matrix analysis that includes the mat- 
ter and baryon density parameters and Qb, fs and 
the non-Gaussian factor G. The results for the l-cr errors 
on G, assuming priors from CMB, Large-Scale Structure 
(LSS) and supernovae (SN) observations, are AG ~ 2 
and AG ~ 0.1 for the Bolocam and Planck experiments 
respective ly. 

Finally, iSadeh et all l|200fiD apply the same extended 



PS formalism to the y'L non-Gaussian model l) White! 
[1999: Kovam a et al.l Il999() . Here, however, much at- 
tention is devoted to highly non-Gaussian models, e.g., 
with m = 1 and 2, which are already excluded by mea- 
surements of the gala xy bispectrum in the PSCz survey 
ijFeldman et alJl200l|) . 

2.2. The non-Gaussian mass function 

In our analysis we will make use of the EPS approach, 
Eq. jni , since it can be more easily implemented (once the 
probability function PNciy) is known) and avoids prob- 
lems with small regions of the parameter space where the 
MVJ expression for the mass function, Eq. © , is beyond 
its limits of validity. For most of the cases considered in 
section|21 however, we performed the analysis using both 
approaches, finding almost identical results. 

Since the PS and EPS expressions are known to differ 
by up to 25% from N-body results, we use the EPS non- 
Gaussian mass function only to model departures from 
the Gaussian case; for the latter we use an analytic mass 
function fit to the N-body results. Specifically, we con- 
sider the non- Gaussian mass function n(z, M, Jnl) to be 
given by the product 

n{z, M, /jvl) = naiz, M) Fng{z, M, /nl), (8) 

where nQ{z,M)^ corresponding to the Gaussian case, is 
the fit to N-body simulations provided bv i.Tenkins et al.l 
(l200l . 



nG{z,M)dM 
exp [ 



d(7M 



0.301 

Mctm dM 



0.64-log[D(z)aM]|'-''] , (9) 

where D{z) is the linear growth factor computed by solv- 
ing the differential equation governing structure evolu- 
tion. The non-Gaussian factor Fmg{z, M, Jnl) is derived 
from the EPS mass function and simply given by 



FNGiz,M, fNL 



_ nEPsiz,M, /nl) 



(10) 



nps{z,M) 

where nps is the Gaussian PS mass function. Note that 
for fpfL = we have ueps = nps- 

It can be easily shown that the predictions of the EPS 
and MVJ methods are very close by comparing them 
for relevant values of the parameter //vl- In Fig. ^ we 
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Fig. 1. — Uncertainty on the mass function n{M,z) due to 
non-Gaussianity expressed as the ratio between the non-Gaussian 
to the Gaussian mass function at four different redshifts (0, 0.5, 
1 and 1.5) for different values of /jvl- The inner continuous lines 
corresponds to the (2-ct) limits —27 < /jvl < 121 derived from 
the W MAP (first year) constraints on /jvl by ICreminelli et all 
120061) while the outer continuous lines corresponds to the lim- 
its —243 < /ivL < 337 from the expected SDSS galaxy bispectrum 
constraints I Scoccimarro ct al. 2004} computed with the EPS ap- 
proach, Eq. J^. The dashed linos, almost coincident with the con- 
tinuous ones, correspond to the same limits computed by means of 
the MVJ formula, Eq. g}. 



plot the ratio of the non-Gaussian mass function to the 
Gaussian one at different redshifts and as a function of 
the mass M for the 2-a limits 

-27 < fNL < 121, 

obtained fro m the bispectrum analysis of the WMAP 1- 
year data bv lCreminelh et aD(j2006') . yielding constraints 
that are slightly tighther than but consist ent with those 
obtained from the WMAP 3-year data bv lSpergel et alJ 
(12006^) . We plot as well the limits 

-243 < fNL < 337, 

corresponding to the l-cr error A/jvl — 145 ex- 
pected from measurements of the galaxy bispectrum 
in the Sloan Digital Sky Survey (SDSS) main sample 
l)Scoccimarro et ;il]l2004(i . Notice that in the latter case 
we are assuming here as fiducial value for the non- 
Gaussian parameter //vi = 47, i.e. the maximum likeli- 
hood value of the cited WMAP analysis. The continu- 
ous lines denote therefore the allowed region computed 
by means of the EPS formula, Eq. lO, while the dashed 
lines corresponds to the same quantity determined as- 
suming the EPS formalism of Eq. ©. 

The non-Gaussian PDF PNGiu) that appears in the 
EPS formula, Eq. is measured from realizations on a 
512'^ grid in a box of 1 h^^ Gpc side of the curvature per- 
turbation $ described by Eq. ^ in terms of the Gaussian 
field generated in Fourier space with a scale invariant 
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power spectrum. The field <i> is then converted into the 
mass density field in Fourier space by means of Poisson 
equation an d a transfer f unction c ompu ted by the CMB- 
FAST code l|Seliak fc Za ldarriagal lT99^ . smoothed on a 
R = 4/i~^Mpc scale and the probability distribution is 
finally measured in position space. We determine a prob- 
ability function for a set of values of Jnl from —500 to 500 
and then interpolated to the desired value. We assume, 
in all cases, the fiducial cosmology described below. In 
order to reliably estimate the tail of the distribution, sev- 
eral thousands realizations were needed for each value of 
Jnl- For high values of the mass M and redshift, corre- 
sponding to extremely rare events (7-ct), the probability 
distributions could not be properly determined, as can be 
seen from the limited plots in the lower panels in Fig. ^ 
None of the results in the paper, however, is sensitive to 
this cut-off. 

It is evident from the figure that the difference be- 
tween the two approaches, for this non-Gaussian model, 
is small, essentially noticeable just for large values of 
Inl- This is due to the relatively mild dependence on the 
smoothing scale R of the reditced skewness s^^r = S^^hUr 
for our non-Gaussian model, as we tested as well by 
choosing different smoothing lengths for the probabil- 
ity distributions measured from the realizations. On the 
other hand, the close results obtained by the different 
methods show how the degree of non-Gaussianity cur- 
rently allowed can be described by the first moments of 
the primordial distribution, if not by the skewness alone. 
It is worth stressing, however, that both prescriptions 
for the non-Gaussian mass function, even when limited 
to modeling deviations from the Gaussian case, need to 
be properly tested against N-body simulation s. New re- 
sults in this direction will soon be available l|Matarresel 

By means of the mentioned measured probability func- 
tions, in the framework of the EPS approach, it is also 
possible to translate the constraints on the parame- 
ter ff4L into constraints on the parameter G defined in 
Eq. Q. As an example, the WMAP 1-cr error A/at^ = 37 
iGrcminclli ct al. 2006) corresponds to AG ~ 0.06 while 
the expected l-tr error A/jv^ — 145 from SDSS galaxy 
bispectrum measurement corresponds to AG ~ 0.25. 

2.3. Likelihood analysis 

In this section we describe the likelihood analysis we 
use to obtain our results. We will consider two simple 
models depending on four and five parameters. In ad- 
dition to //VL we consider the matter density parame- 
ter rim, fluctuation amplitude parameter erg and we will 
separately consider the cases of dark energy with either 
a constant equation of state parameter [w) or a time- 
varying equation of state described by two parameters 
(wo and Wa)- In all cases we assume a spatially flat cos- 
mological model for simplicity. 

The fiducial values assumed for the likelihood analysis 
are given in Table ^ Since the expected number of ob- 
servable clusters is highly dependent on the value of erg, , 
for the four parameter model we perform the analysis as- 
suming as well the lower value CTg = 0.75, while in every 
other case we assume tis = 0.9. The choice of the fidu- 
cial value /jvL = 47 for the non-Gaussian parameter does 
not substantially affect any of the results of the present 
work. 



TABLE 1 

Fiducial values for the cosmological and non-Gaussian 
parameters. 



Parameter 




Fiducial value 




matter density 


0.27 


0"8 


galaxy fluctuation amplitude 


0.9 (0.75) 




dark energy equation of state 


-1 


Wa 


dark energy equation of state 





Jnl 


non-Gaussian parameter 


47 




scalar spectral index 


1 


h 


Hubble parameter 


0.72 




physical baryon density 


0.0232 




dark energy density 


1 - nm 



Unless otherwise stated, we consider an ideal survey 
with limiting mass M/™ = 1.75 X 10^'' h-^ Mq and with 
a sky coverage of 4000 deg^ {fsky — 10%) out to a maxi- 
mum cluster redshift of 1.5, corresponding to the expec- 
tations for the SPT and DES projects. For our fiducial 
model with erg — 0.9 and /jv^ = 47, this yields a to- 
tal of 21,000 clusters in 15 redshift bins; if we had in- 
stead chosen Jnl = for the fiducial model, we would 
obta in 20 , 000 clusters, consistent with earlier estimates 
tWang et al. 2004). 

We study the dependence on cosmology and on the 
constant Jnl of the total number and mass distribution 
of clusters above a certain fixed, i.e. redshift independent, 
threshold mass Mum and explore the degeneracies intro- 
duced by varying non-Gaussian initial conditions. While 
the redshift dependence of the threshold mass should be 
included when making precise predictions for a given sur- 
vey, this dependence is weak for SZ-selected cluster sam- 
ples; as a result, our neglect of such dependence here will 
not significantly affect our conclusions. 

The total number of clusters with mass M above Mum, 
per unit redshift, is given by 



dN 
dz 



dzdfl 



n{z,MjNL)dM (11) 



where 



dV 
dzdQ 



H{z) 



dz' 
H{z') 



(12) 



is the cosmology-dependent volume factor for flat models 
and Afi is the solid angle subtended by the survey area. 

We show in Fig.Elthe sensitivity of the total number of 
clusters above Mum per unit redshift and unit area (up- 
per panels) and of the comoving number density (lower 
panels) to different values of the non-Gaussian parame- 
ter //VI compared with the sensitivity to different values 
of the dark energy equation of state parameter wq and 
of the fluctuation amplitude parameter erg . 

Fig. El (upper left panel) shows that varying fj^^ over 
the range allowed by current CMB observations yields 
changes in the cluster counts comparable to a 10% vari- 
ation in the dark energy equation of state parameter w. 
However, the upper right panel of Fig. El shows that the 
redshift dependence of the mass function variations due 
to non-Gaussianity are different from the variations due 
to changes in w. This is essentially due to the fact that 
w affects both the mass function and the volume factor. 
On the other hand, the redshift dependence of variations 



5 





o 2 



.1 11 1 1 1 
: z=0 


1 ' 'I 




//" 




// : 




// : 




^ / : 




1 . r 




10" 10'= 
M [h-'Mg 



10'* 10" 
M [h-'M^] 



Fig. 2. — Cluster counts per unit redshift (upper panels) and 
comoving cluster density (lower panels) as a function of redshift for 
different values of the non-Gaussian parameter /jvi (/jvl = TlOO, 
continuous lines), of the dark energy equation of state parameter 
w {w = —1.1 and —0.9, short-dashed lines) and of as (erg = 0.85 
and 0.95, long-dashed lines) compared to the fiducial case (dotted 
line) with erg = 0.9, ui = —1 and /jvi = 0. Assumes the mass limit 
Miirr, = 1.75 X 1014 h-^MQ. 



Fig. 3. — The mass function ■n{M, z) as a function of the mass M 
at z = (upper panels) and z = I (lower panels) for different values 
of the non-Gaussian parameter /]vl (continuous lines), of the dark 
energy equation of state parameter wq (short-dashed lines) and of 
(Tg (long-dashed lines) compared to the fiducial case (dotted line) 
with (Tg = 0.9, w = ~l and /ati, = 0. 



due to changes in //vi appears more similar to variations 
induced by changes in ag , so we expect a stronger degen- 
eracy between these two parameters. 

In Fig. 13 we show the sensitivity of the mass function 
n[M, z) on the same parameters, this time as a func- 
tion of the mass M for z = (upper panels) and z = 1 
(lower panels). In this case the behavior of the cluster 
density as we vary //vl and cts is quite different. One can 
clearly see how non-Gaussianity is particularly significant 
for the high-mass tail of the distribution. This fact sug- 
gests that it might be relevant to consider a likelihood 
analysis that takes into account the full functional shape 
of the mass function by dividing the observable clus- 
ters in mass bins (see, fo r instance, Hu 2003; Lima & Hu 
arian fc Bernstein 20 06). In this way one might 
expect to mitigate the /nl-^s degeneracy evident from 
Fig. [3 and better study the possibility of constraining 
non-Gaussianity with cluster surveys. 

In the next section we will consider the two cases of an 
analysis involving a single mass bin defined by M > Mum 
and of an analysis with several mass bins. The likelihood 
function is based on the assumption of Poisson statistics 
for the cluster number measurements in each redshift bin, 
so that, for the single mass bin case we have 



ln£ 



[N, In N*-N*- N, In N, + N,] (13) 



where N^^^ is the total number of redshift bins, N* is 
the fiducial number count in the i-th redshift bin, Ni is 
the number count of clusters in the i-th redshift bin for 



the specific model. Throughout the paper, we consider 
iVjQj = 15 redshift bins with width Az = 0.1 out to 
z = 1.5. 

In the case of multiple mass bins, the likelihood func- 
tion is of the form 



i=i j=i 

-N,j\nN,j+N,,] (14) 



with Ni^f — 10 being the total number of mass 
logarithmically spaced from to 



bins 

X 

10^^ /i""'^ Mq; here TVy is the number of model clusters 
in the i-th redshift bin and j-th mass bin. 

The marginalization of the likelihood functions is per- 
formed on a regular grid with a varying number of points 
chosen to optimize the sampling of the parameter space. 

We do not include systematic errors which a real sur- 
vey will encounter, including uncertainties in the cluster 
mass - observable relation (e. g. Scliak 2002; Lcvinc et al] 
2002t iPierpaoli et alJ 120031: Maiumdar fc Mohr .2005 
2003 lHn2003t [Lima fc HiJ 120041 iFrancis et all 1200^ 

" " " :ooa ~ " 



Lima & Hul 1200,4 Kravtsov et aimOOfil) and in cluster 
redshift determination fe.g. lHuterer et al.ll2004(l . We ex- 
cluded as wel l statistical unc ertainties related to sample 
variance (e.g. lHu fc Kravts'ov^2003) or theoretical uncer- 
tainties in the clust er mass function, and its cosmologi- 
cal depend ence (e .g. *Heitmanji_jt_ajj2005t IWarren et all 
12006: Ree d et al.l l2006: Cro cce et al.fl2006l) . Finally, we 
note that other degeneracies with para meters such as 
those describing spatial curvature (| Abbott et al.ll2005|) 
and the effect of massive neutrinos on the dark matter 
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power spectrum l|Huterer fc Linderll2006(l might be rele- 
vant for future high precision analyses. 

3. RESULTS 

In this section we estimate the impact of marginalizing 
over the non-Gaussian parameter /^tl on the determina- 
tion of the dark energy equation of state as well as on 
two other relevant cosmological parameters such as the 
matter density flm and fluctuation amplitude as- We 
will separately consider the case of a dark energy equa- 
tion of state determined by a single parameter (w) and 
the case of a two-parameter description of a time- varying 
equation of state {wq and Wa)- 

We derive the marginalized errors on the parameters 
with fixed /nl — 47 (no marginalization) and with three 
different Gaussian priors on /jvl, two corresponding to 
the constraints from CMB bispectru m measurements ex- 
pected from the Planck experiment ijKomatsu & SoergeJ 
UpOl; Liguori et al. 2006) an d mea sured in the WMAP 
experiment IjCreminelli et 311120061) with 



fNL =47 ±5 (1-CT, Planck) 



(15) 



and 



/jvl = 47 ±37 (l-fT, WMAP), (16) 

and a third corresponding to the expected constraints 
from the analysis of the SP SS main sample galaxy bis- 
pectrum ( Scoccimarro et al.i.200 4). 

/ArL=47±145 (1-cr, SDSS forecast). (17) 

This last case is motivated by a possible strong scale- 
dependence of primordial non-Gaussianity, not captured 
by the model defined by Eq. that could result in a 
stronger non-Gaussian effect at smaller scales, thereby 
escaping the CMB constraints. As a rough estimate 
of the smallest scale probed by the mentioned experi- 
ments, we notice that for WMAP, the maximum mul- 
tipole Imax — 1000 corresponds to ~ 50/i~^Mpc while 
Planck is expected to probe a scale three times smaller; 
in the SDSS case a maximum comoving wavenumber 
k,nax — 0.3/iMpc~^ corresponds to 20/i~^Mpc. The 
typical scale probed by clusters is about 5 to 10 h~-^ Mpc, 
with the most massive clusters approaching the lowest 
scale probed by Planck. 

In all the different cases considered we include as well 
the results obtained with two, independent, Gaussian pri- 
ors on Qrn and cg with errors roughly corresponding to 
the knowledge provided by WMA P observations for a 
ACDM model l|Spergel et alJl2006j) in combination with 
other probes, such as, for example, the LSS power spec- 
trum, 

<T8 = 0.9 ±0.05, and a™ = 0.27 ± 0.035, (18) 

and by future con straints from Pl anck in combination 
with other probes l|PL ANCKll2"n(i^ . 

cTg = 0.9 ±0.01, and = 0.27 ± 0.0035. (19) 

As an extreme example, in the last two lines of tables, we 
give results corresponding to fixing flm and as , studying 
therefore a likelihood function for the dark energy pa- 
rameters and Jnl alone. 

We caution that these priors have a purely illustrative 
significance and are chosen here for the sake of simplic- 
ity. A proper treatment of external data sets, which is 



TABLE 2 

Expected cosmological errors (1 — o-) from the cluster 

SURVEY FOR THE 4-PARAMETER (Qrn, erg, W, /jvl) ANALYSIS. THE 
PERCENTAGES IN PARENTHESES EXPRESS THE INCREASE IN THE 
ERROR WITH RESPECT TO THE CASE WITHOUT MARGINALIZATION 
ON /jvl (A/jvl = 0). We assume a fiducial as = 0.9 AND ONE 
MASS BIN DEFINED BY M > Murn = 1-75 X lO^** Mq. 



prior : 


A fNL = 





A/nl = 5 


A/nl = 37 


A/jvt = 145 


No priors 


on 


and 


CTS 






A-u; 




0.045 




0.045 (0%) 


0.049 (9%) 


0.079 (76%) 


Ar2„, 




0.0085 




0.0085 (0%) 


0.0085 (0%) 


0.0087 (2%) 


Aa-s 




0.0051 




0.0052 (2%) 


0.0083 (63%) 




AfNL 








5.0 


37 


123 


Gaussian 


priors; 




= 0.27 ± 0.035, 


CTs = 0.9 ± 0.05 




A-u; 




0.044 




0.044 (0%) 


0.048 (9%) 


0.076 (73%) 


Ar2„ 




0.0082 




0.0082 (0%) 


0.0082 (0%) 


0.0083 (1%) 


Acts 




0.0050 




0.0050 (0%) 


0.0081 (62%) 


0.0205 (310%) 


A/jvt 








5.0 


36 


113 


Gaussian 


priors; 




= 0.27 ± 0.0035, CTg = 0.9 ± 0.01 




Aw 




0.023 




0.024 (4%) 


0.031 (35%) 


0.042 (83%) 


An,„ 




0.0032 




0.0032 (0%) 


0.0032 (0%) 


0.0032 (0%) 


Acts 




0.0021 




0.0023 (10%) 


0.0055 (160%) 


0.0091 (330%) 


A/JVL 








5.0 


31 


54 


Fixed Sl„ 


= 0.27 and 


CTs = 0.9 






Aw 




0.0172 




0.0177 (3%) 


0.0184 (9%) 


0.0184 (9%) 


A/jvt 








3.8 


5.6 


5.7 



beyond the scope of this paper, would naturally involve 
the parameters covariance and it would affect directly the 
dark energy parameters as well. On the other hand, even 
rigorous analyses of CMB or LSS galaxy power spectra 
would probably be insensitive to the non-Gaussian pa- 
rameter /jvl. 

3.1. 1-parameter Dark Energy equation oj state 

The main results in this paper are shown in Table 
where we present the expected 1-a errors from the cluster 
survey for the three parameters fim, o'8 and w with no 
marginalization on /jvl (A/jvl = 0) and with a marginal- 
ization that includes the three Gaussian priors discussed 
above (A/jvl = 5, 37 and 145), assuming in this case a 
fiducial as — 0.9. The percentages in parentheses express 
the increase in the error with respect to the case without 
marginalization on /jvl • Although the derived marginal- 
ized likelihood for single parameters are quite close to 
Gaussian functions, the estimated errors reported in Ta- 
ble 12 as in the following ones, are given, for clarity, by 
the mean between upper and lower errors. 

The major conclusion from Table |2] is that inclusion of 
a possible non-Gaussian component at the level allowed 
by present and future CMB constraints will not have an 
appreciable impact on the determination of dark energy 
parameters from cluster surveys. On the other hand, us- 
ing only a single cluster mass bin (i.e., no information 
about the shape of the cluster mass function) and the 
WMAP prior on /jvl, the inclusion of non-Gaussianity 
degrades the cluster constraint on as by 50%. Moreover, 
using only the projected SDSS bispectrum constraint on 
//VL, we do see degeneracies between /jvl and dark en- 
ergy: the error on w from clusters increases by ~ 70%, 
and the error on as grows by a factor of more than three 
compared to the purely Gaussian case. In all cases, the 
determination of stays largely unaffected. 



7 



0.95 - 



oo 

b 0.9 



0.85 
-0.8 

-1.2 

250 


-250 



0.25 



fi^i^=47 fixed 



WMAP prior: fNL=4'J' ± 37 

_ _ LSS prior: fNi,=47 ± 145 






It. ■ I 



I ■ .t ■ I 



all contours 95% CL 



- 0.8 



0.6 
d 0.4 



n. 



0.3 0.85 0.9 0.95 
''a 



-1.2 



-0.8 -250 



w 



250 



NL 



Fig. 4. — Forecast marginalized likelihoods and 95% C.L. contour plots from the cluster survey for the 4-parameter (Qm, eg, w and 
Jnl) model. In each contour plot, the other two parameters are marginalized over; in each likelihood plot, the other three parameters are 
marginalized over. We assume a fiducial ACDM model with erg = 0.9 and use one mass bin defined by M > = 1.75 X 10^* Mq. 



The expected degeneracy between /jvl and erg is ev- 
ident from the two-parameter 95% C.L. contour plots 
shown in Fig. 01 where we marginaUzed over the re- 
maining parameters. The same overall behavior is ob- 
served when we impose priors on and erg ■ In Fig- 
the marginalized likelihood plot for w shows that inclu- 
sion of the WMAP prior on //vl leads to essentially the 
same dark energy sensitivity for the cluster survey as 
one would have by fixing /at^ (i.e., by not including non- 
Gaussianity) . 

As far as the constraints on non-Gaussianity are con- 



cerned, the information provided by the cluster likelihood 
on /jvi adds little to that from the CMB, while it mildly 
improves upon the non-Gaussianity constraints expected 
from the galaxy bispectrum. Even in the ideal case of 
perfect knowledge of Qm and erg, the expected error on 
Jnl is ^/nl — 6, of the same order of the expected error 
from Planck. 

We performed the same analysis for a fiducial erg = 
0.75, since this lower value has been recently suggested 
by C^'IB (SDcrgel ct al. 20 06il and cluster observations 
iGladders et a,lJl2nn4 iDahfFn2nnfil^ : as is well known, a 
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TABLE 3 

Same as Table |^ but with fiducial erg =0.75. 



prior: 


A/jvt = 


AfNL = 5 


AfNL = 37 


AfNL = 145 


No priors 


Aw 




0.079 


0.079 (0%) 


0.083 (5%) 


0.124 (57%) 


An„ 




0.0140 


0.0140 (0%) 


0.0140 (0%) 


0.0144 (3%) 


Acts 




0.0076 


0.0076 (0%) 


0.0100 (32%) 


0.0238 (210%) 


A/nl 




- 


5.0 


37 


128 


Gaussian 


priors; Qj^ 


= 0.27 ± 0.035, 


as = 0.75 ± 0.05 








0.073 


0.073 (0%) 


0.079 (8%) 


0.119 (63%) 


An„ 




0.0129 


0.0129 (0%) 


0.0129 (0%) 


0.0129 (0%) 


Ao-g 




0.0070 


0.0070 (0%) 


0.0094 (34%) 


0.0212 (200%) 








5.0 


37 


118 


Gaussian 


priors; Qrn 


= 0.27 ± 0.0035, as = 0.75 ± 0.01 


Aw 




0.036 


0.036 (0%) 


0.046 (28%) 


0.063 (75%) 


An„ 




0.0033 


0.0033 (0%) 


0.0034 (3%) 


0.0034 (3%) 


Ao-g 




0.0021 


0.0023 (10%) 


0.0054 (160%) 


0.0091 (330%) 


A/jvi 






5.0 


32 


56 


Fixed f2„ 


= 0.27 and 


as — 0.75 






A-u; 




0.030 


0.031 (3%) 


0.034 (10%) 


0.034 (10%) 


A/ivt 






4.2 


7.5 


7.6 



lower clustering amplitude reduces the number of ex- 
pected clusters and thereby reduces the constraining 
power of cluster surveys. The results are given in Ta- 
ble 121 In this case the total number of clusters for the 
fiducial model is about 6, 000; as a consequence, the cos- 
mological constraints from the cluster survey are weaker 
than for the high-erg model. However, the relative impact 
of the marginalization over primordial non-Gaussianity 
is reduced. This result can be expected since the effect 
of imposing the same priors on //vl is relatively smaller 
when the cosmological errors on the other parameters for 
the fixed /nl case increase. 

To further illustrate the dependence of the results on 
survey parameters, in Table 0] we show the constraints 
obtained when the threshold cluster mass is reduced to 
Mii^ = 1 X 10^'*/i"^Mo. This lower threshold may 
be achieved, e.g., by supplementing SZ cluster detec- 
tion with optical cluster selection using the red galax y 
sequence fe.g- .Gladders et al.ll2006t iKoester et al.ll2006j) . 
In this case, the 4,000 deg^ survey to z = 1.5 includes 
about 75, 000 clusters, and the forecast cosmological pa- 
rameter errors (without non-Gaussianity) are smaller by 
almost a factor of two than for the case with larger Mum 
considered above. For this more sensitive cluster sur- 
vey, the impact on cosmological parameters of marginal- 
izing over /jvL is correspondingly larger: while the impact 
on dark energy remains small, including non-Gaussianity 
with the WMAP prior expands the error on erg by more 
than 100%. 

As already discussed in the previous section, the degen- 
eracy between ag and Jnl could be partially reduced by 
introducing a number of cluster mass bins and using the 
information contained in the shape of the mass function. 
In Table |S1 we present the results for an analysis with a 
fiducial erg = 0.9 and M^m = 1.75 x lO^^h^^MQ as in 
Table 121 but subdividing the clusters into ten mass bins 
and using the likelihood function defined in Eq. (|14|) . As 
the last column in Table [S] indicates the main effect of 
including mass bins is that the cluster constraint on the 
non-Gaussian parameter /jvl becomes stronger than that 



TABLE 4 

Same as Table^but with Murn = 1 x 10^'' Mq. 



prior ; 


AfNL = 


AfNL = 5 


AfNL = 37 


AfNL = 145 


No priors 


on Sl^ and 








A-u; 




0.026 


0.026 (0%) 


0.030 (15%) 


0.052 (100%) 


An,„ 




0.0050 


0.0050 (0%) 


0.0050 (0%) 


0.0052 (4%) 


Acts 




0.0031 


0.0032 (3%) 


0.0066 (110%) 


0.0186 (500%) 


A/nl 




- 


5.0 


36 


113 


Gaussian 


priors; 


= 0.27 ± 0.035, 


as = 0.9 ± 0.05 




Aw 




0.026 


0.026 (0%) 


0.030 (15%) 


0.050 (92%) 


An,„ 




0.0050 


0.0050 (0%) 


0.0050 (0%) 


0.0051 (2%) 


Aas 




0.0030 


0.0031 (3%) 


0.0066 (120%) 


0.0174 (480%) 


AfNL 






5.0 


36 


106 


Gaussian 


priors; 


= 0.27 ± 0.0035, as = 0.9 ± 0.01 




Au; 




0.017 


0.017 (0%) 


0.023 (35% 


0.032 (88%) 


An,„ 




0.0028 


0.0028 (0%) 


0.0028 (0%) 


0.0028 (0%) 


Acts 




0.0018 


0.0020 (11%) 


0.0051 (180%) 


0.0087 (380%) 


A/nl 






5.0 


32 


56 


Fixed n,„ 


= 0.27 and 


as = 0.9 






Aw 




0.0100 


0.0102 (2%) 


0.0102 (2%) 


0.0102 (2%) 


AfNL 






3.1 


4.0 


4.0 



from the SDSS galaxy bispectrum. Even without com- 
bining with external data sets one can reach a 1 — cr error 
of A/tvl ~ 50, not too far from current limits from CMB 
observations. Further study would be needed to deter- 
mine if this conclusion remains when realistic uncertain- 
ties in the cluster mass-observable relation are included 
in the analysis. 

3.2. 2-parameter Dark Energy equation of state 

Finally we consider the case of a time-vary ing dark 
energy equa tion of state l|Chevallier &: Polarski 200ll 
ILindeHl200^ . 

w{a) = Wo + (l - a)wa, (20) 

adding the parameter Wa to the likelihood analysis stud- 
ied so far. In this case the strong degeneracy between 
Wo and Wa enlarges considerably the region of parameter 
space that has to be covered for the likelihood function 
evaluation, including unphysical regions where the com- 
bination wq + Wa, representing the equation of state at 
large redshift, takes large positive values. To avoid such 
cases we impose, by hand, a Gaussian prior on the value 
of f2„j(z), requiring in particular l — Q^iz) < 0.01 at 1-cr 
at z = 30, the initial redshift considered for the numer- 
ical solution to the differential equation governing the 
growth factor D{z). This ensures that the Universe is 
matter-dominated at early times as required by struc- 
ture growth. 

In Table we present the derived 1-cr errors on the 
five parameters fim, erg, wq and Wa with and without 
marginalization on /nl and assuming a single mass bin 
defined by M > Mu^ = 1.75x10" h'^ Mq. Since in this 
case the uncertainties on the parameters, which are sensi- 
tive to the strong WQ-Wa degeneracy, are much larger than 
in the previous case, the effect of the marginalization on 
Inl with a CMB prior is even smaller than in the case 
of time-independent w. As before, however, marginal- 
ization over f^L with only the galaxy bispectrum prior 
substantially increases the error on erg. 

As a final example, in Table we consider parameter 
constraints in the time-varying dark energy cosmology 
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using the ten cluster mass bins. As noticed earlier in the 
case of the 4-paranieter analysis, the ag-fNL degeneracy 
is significantly reduced and the expected constraints on 
non-Gaussianity are still of the order of A/jvi ~ 50. 

4. CONCLUSIONS 

The success of the ACDM standard cosmological model 
in recent years has been nothing short of spectacular. 
Upcoming surveys will either continue to confirm this 
model and constrain its parameters with unprecedented 
accuracy, or they will uncover discrepancies which will 
point the way toward improvements in our understand- 
ing of fundamental physics. Two questions addressing 
cosmology beyond the standard model that have been 
the subject of substantial attention in recent years are: 
what is the nature of the dark energy which is driving the 
accelerated expansion of the universe? and second, are 
fluctuations in the primordial matter distribution Gaus- 
sian, and therefore consistent with the predictions of the 
simplest inflationary models? In this paper, we obtain 
a rough estimate of the success that one of the most 
promising cosmological probes, galaxy cluster counts, is 
likely to have in answering these fundamental questions. 

We have assumed an ideal cluster survey with survey 
parameters expected for the upcoming SPT and DES 
projects. Our fiducial cosmological model includes both 
dark energy and primordial non-Gaussianity using popu- 
lar parameterized models, w and Wa for the former, and 
Inl for the latter. Cluster number counts as a func- 
tion of mass, redshift, and cosmolo gy were estimated us - 
ing a standard fit to simulations l|Jenkins et alJ 120011) . 
which we adjusted to allow for mild non-Gaussian ini- 
tial conditions, and all clusters above a threshold mass 
were considered to be "found" by our fiducial survey. 
We then performed a simple likelihood analysis on the 
cluster counts using priors from current WMAP and ex- 
pected Planck and SDSS constraints on non-Gaussianity 
as well as approximate priors on the two other relevant 
cosmological parameters from other present and future 
data sets. 



TABLE 5 

Same as Table|^but using 10 cluster mass bins. 



prior : 


A/jvL = 


AfNL = 5 


AfNL = 37 


AfNL = 145 


No priors 


on and 


CT8 






Aw 




0.044 


0.044 (0%) 


0.046 (5%) 


0.048 (9%) 


An„ 




0.0082 


0.0082 (0%) 


0.0083 (1%) 


0.0085 (4%) 






0.0049 


0.0050 (2%) 


0.0077 (57%) 


0.0103 (110%) 








5.0 


29 


45 


Gaussian 


priors; 


= 0.27 ± 0.035, 


CTg = 0.9 ± 0.05 




A-u; 




0.043 


0.043 (0%) 


0.045 (5%) 


0.048 (12%) 


An„ 




0.0079 


0.0079 (0%) 


0.0080 (1%) 


0.0082 (4%) 


Acts 




0.0048 


0.0049 (2%) 


0.0075 (56%) 


0.0100 (110%) 


A/nl 






5.0 


29 


45 


Gaussian 


priors; 0„i 


= 0.27 ± 0.0035, CTs = 0.9 ± 0.01 




Aw 




0.023 


0.023 (0%) 


0.029 (26%) 


0.034 (48%) 


An„ 




0.0032 


0.0032 (0%) 


0.0032 (0%) 


0.0032 (0%) 


Acts 




0.0021 


0.0022 (5%) 


0.0048 (130%) 


0.0063 (200%) 


A/jvt 






5.0 


26 


35 


Fixed 0„ 


= 0.27 and 


CTS = 0.9 






A-u; 




0.0166 


0.0173 (4%) 


0.0181 (9%) 


0.0181 (9%) 








3.7 


5.5 


5.5 



TABLE 6 

Expected cosmological errors (la) from the cluster 

SURVEY FOR THE 5-PARAMETER (Qm, erg, Wq, Wa, Jnl) ANALYSIS. 
We ASSUME A FIDUCIAL CTg = 0.9 AND ONE MASS BIN DEFINED BY 

Mu^ = 1.75 X lOi4/i-iM0. 



prior; 


AfNL = 





A/ni, = 5 


AfNL = 37 


AfNL = 145 


No priors on 


and 


CTs- 






Awo 


0.195 




0.195 (0%) 


0.196 (1%) 


0.208 (7%) 


Awa 


0.73 




0.73 (0%) 


0.74 (1%) 


0.87 (19%) 


An^ 


0.0156 




0.0156 (0%) 


0.0158 (1%) 


0.0185 (19%) 


Acts 


0.0086 




0.0087 (1%) 


0.0114 (32%) 


0.0296 (240%) 


A/nl 






5.0 


37 


137 


Gaussian priors; 




= 0.27 ± 0.035, 


CTg = 0.9 ± 0.05 




Awo 


0.178 




0.178 (0%) 


0.178 (0%) 


0.183 (3%) 


AWa 


0.68 




0.68 (0%) 


0.68 (0%) 


0.73 (7%) 


An^ 


0.0140 




0.0140 (0%) 


0.0141 (1%) 


0.0151 (8%) 


Acts 


0.0078 




0.0078 (0%) 


0.0105 (35%) 


0.0240 (210%) 


A/nl 






5.0 


37 


118 


Gaussian priors; 




= 0.27 ± 0.0035, CTs = 0.9 ± 0.01 




Awo 


0.082 




0.082 (0%) 


0.85 (4%) 


0.89 (9%) 


AlUa 


0.43 




0.43 (0%) 


0.43 (0%) 


0.43 (0%) 


Ari„, 


0.0034 




0.0034 (0%) 


0.0034 (0%) 


0.034 (0%) 


Acts 


0.0021 




0.0023 (10%) 


0.0055 (160%) 


0.0091 (330%) 


A/nl 






5.0 


31 


54 


Fixed n,„ = 0.27 and 


CTs = 0.9 






Awo 


0.068 




0.069 (1%) 


0.071 (4%) 


0.071 (4%) 


Awa 


0.37 




0.039 (5%) 


0.40 (8%) 


0.40 (8%) 


A/lVL 






3.9 


6.1 


6.2 



Our principal conclusion is that dark energy con- 
straints are in all cases not substantially degraded by 
primordial non-Gaussianity when the model parameter- 
ized by the constant //vi and current limits from CMB 
observations are assumed. This is true despite the fact 
that variations in f^L close to current uncertainties in- 
duce differences in the mass function comparable in mag- 
nitude to variations of 10% in the dark energy parameter 
w. A stronger degeneracy is observed instead between 
fNL and fjg; in this case, the expected errors on erg from 
future cluster surveys can be noticeably affected when 
non-Gaussianity is included in the analysis. 

A secondary conclusion is that the cluster survey it- 
self might have sufficient statistical power to provide a 
valuable cross check on any detection or non-detection of 
primordial non-Gaussianity in CMB experiments, partic- 
ularly when information on the cluster distribution as a 
function of the mass is taken into account. 

However, we must emphasize that we have not at- 
tempted to include in our analysis any of the system- 
atic and statistical errors in the clusters mass determi- 
nation, which are likely to cause trouble for real surveys, 
as well as uncertainties on the predictions for the mass 
function, and our results must be interpreted with this 
in mind. We believe that our principal result should be 
quite robust, since any significant increase in the error 
budget will reduce constraining power on dark energy 
parameters and de-emphasize the confusion caused by 
any non-Gaussian initial conditions. On the other hand, 
the effectiveness of clusters as a cross check of primordial 
non-Gaussianity estimates from the CMB could be dra- 
matically worsened, and should therefore be the subject 
of future work. 
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TABLE 7 

Same as TableIHIbut using 10 cluster mass bins. 



prior: 


A/jvL = 


A/nl = 5 


A/jvi, = 37 




No priors on 0„i and 


o"8. 






Awq 


0.180 


0.180 (0%) 


0.186 (3%) 


0.197 (9%) 




0.68 


0.68 (0%) 


0.71 (4%) 


0.78 (15%) 




0.0141 


0.0141 (0%) 


0.0151 (7%) 


U.Ul/U \Z1 /(I ) 


Ao-g 


0.0078 


0.0079 (1%) 


0.0113 (45%) 




A/jvi 




5.0 


31 


51 


Gaussian priors; 


= 0.27 ± 0.035, as = 0.9 ± 0.05 




AlDo 


0.167 


0.168 (1%) 


0.170 (2%) 


0.175 (5%) 




0.64 


0.64 (0%) 


0.66 (3%) 


0.70 (9%) 


AO^ 


0.0129 


0.0130 (1%) 


0.0136 (5%) 


U.Ul^D t^io/oj 


Ao-g 


0.0072 


0.0073 (1%) 


0.0103 (43%) 




A/jvi, 




5.0 


30 


48 


Gaussian priors; 0„i 


= 0.27 ± 0.003 


5, erg = 0.9 ± 0.01 




Atuo 


0.081 


0.081 (0%) 


0.83 (2%) 


0.84 (4%) 


Awa 


0.42 


0.42 (0%) 


0.42 (0%) 


0.42 (0%) 


An^ 


0.0033 


0.0033 (0%) 


0.0033 (0%) 


0.033 (0%) 


AcTg 


0.0021 


0.0023 (10%) 


0.0049 (130%) 


0.0063 (200%) 


A/nl 




5.0 


26 


35 


Fixed 0„ = 0.27 and 


erg = 0.9 






Awo 


0.067 


0.068 (1%) 


0.070 (4%) 


0.070 (4%) 




0.37 


0.038 (3%) 


0.40 (8%) 


0.40 (8%) 






3.8 


5.9 


6.0 



We thank Eiichiro Komatsu and Eric Linder for helpful 
comments on a earlier draft of the paper. C. V. would 
like to thank Martin White for useful discussions. This 
work was supported by the US Department of Energy at 
the University of Chicago and at Fermilab, by the Kavli 
Institute for Cosmological Physics at the University of 
Chicago and by NASA grant NAG5-10842 at Fermilab. 



REFERENCES 



Abbott, T., et al. 2005, astro-ph/0510346 

Acquaviva, V., Bartolo, N., Matarrese, S., & Riotto, A. 2003, Nucl. 

Phys., B667, 119 
Alishahiha, M., Silverstein, E., & Tong, D. 2004, Phys. Rev.. D70, 

123505 

Amara, A., & Refregier, A. 2004, Mon. Not. Roy. Astron. Soc, 351, 
375 

Arkani-Hamed, N., Creminelli, P., Mukohyama, S., & Zaldarriaga, 

M. 2004, JCAP, 0404, 001 
Bartolo, N., Komatsu, E., Matarrese, S., & Riotto, A. 2004, Phys. 

Rept., 402, 103 
Battye, R. A., & Weller, J. 2003, Phys. Rev., D68, 083506 
Benson, A. J., Reichardt, C, & Kamionkowski, M. 2002, Mon. Not. 

Roy Astron. Soc., 331, 71 
Chen, C, & Szapudi, I. 2006, astro-ph/0606394 
Chevallier, M., & Polarski, D. 2001, Int. J. Mod. Phys., DIO, 213 
Chiu, W. A., Ostriker, J. P., & Strauss, M. A. 1998, Astrophys. J., 

494, 479 

Colafrancesco, S., Lucchin, F., & Matarrese, S. 1989, Astrophys. 
J., 345, 3 

Creminelli, P. 2003, JCAP, 0310, 003 

Creminelli, P., Nicolis, A., Senatore, L., Tegmark, M., & 

Zaldarriaga, M. 2006, JCAP, 0605, 004 
Crocce, M., Pueblas, S., & Scoccimarro, R. 2006, submitted to 

Mon. Not. Roy. Astron. Soc, astro-ph/0606505 
Dahle, H. 2006, submitted to Astrophys. J., astro-ph/0608480 
Dvali, G., Gruzinov, A., & Zaldarriaga, M. 2004, Phys. Rev., D69, 

023505 

Ebeling, H., et al. 1996, Mon. Not. Roy. Astron. Soc., 281, 799 
Feldman, H. A., Frieman, J. A., Fry, J. N., & Scoccimarro, R. 2001, 

Phys. Rev. Lett., 86, 1434 
Francis, M. R., Bean, R., & Kosowsky, A. 2005, JCAP, 0512, 001 
Gladders, M. D., et al. 2006, astro-ph/0603588 
Haiman, Z., Mohr, J. J., & Holder, G. P. 2000, Astrophys. J., 553, 

545 

Heitmann, K., Ricker, P. M., Warren, M. S., & Habib, S. 2005, 

Astrophys. J. Suppl., 160, 28 
Henry, J. P. 1997, Astrophys. J., 489, LI 
Henry, J. P., & Arnaud, K. A. 1991, Astrophys. J., 372, 410 
Holder, G., Haiman, Z., & Mohr, J. 2001, Astrophys. J., 560, Llll 
Hu, W. 2003, Phys. Rev., D67, 081304 
Hu, W., & Kravtsov, A. V. 2003, Astrophys. J., 584, 702 
Huterer, D., Kim, A., Krauss, L. M., & Broderick, T. 2004, 

Astrophys. J., 6l5, 595 
Huterer, D., & Linder, E. V. 2006, astro-ph/0608681 
Jenkins, A., et al. 2001, Mon. Not. Roy. Astron. Soc, 321, 372 
Koester, B., et al. 2006, submitted to Astrophys J. 
Kolb, E. W., Riotto, A., & Vallinotto, A. 2006, Phys. Rev., D73, 

023522 

Komatsu, E., & Spergel, D. N. 2001, Phys. Rev., D63, 063002 
Komatsu, E., et al. 2003, Astrophys. J. Suppl., 148, 119 



Koyama, K., Soda, J., & Taruya, A. 1999, Mon. Not. Roy. Astron. 
Soc, 310, 1111 

Kravtsov, A. V., Vikhlinin, A., & Nagai, D. 2006, submitted to 

Astrophys. J., astro-ph/0603205 
Levine, E. S., Schulz, A. E., & White, M. J. 2002, Astrophys. J., 

577, 569 

Liguori, M., Hansen, F. K., Komatsu, E., Matarrese, S., & Riotto, 

A. 2006, Phys. Rev., D73, 043505 
Lima, M., & Hu, W. 2004, Phys. Rev., D70, 043504 
— . 2005, Phys. Rev., D72, 043006 
Linder, E. V. 2003, Phys. Rev. Lett., 90, 091301 
Lucchin, F., & Matarrese, S. 1988, Astrophys. J., 330, 535 
Lyth, D. H., Ungarelli, C, & Wands, D. 2003, Phys. Rev., D67, 

023503 

Majumdar, S., & Mohr, J. J. 2003, Astrophys. J., 585, 603 
— . 2004, Astrophys. J., 613, 41 
Maldacena, J. M. 2003, JHEP, 05, 013 

Marian, L., & Bernstein, G. M. 2006, Phys. Rev., D73, 123525 
Matarrese, S. 2006, private communication 

Matarrese, S., Verde, L., & Jimenez, R. 2000, Astrophys. J., 541, 
10 

Mathis, H., Diego, J. M., & Silk, J. 2004, Mon. Not. Roy. Astron. 
Soc, 353, 681 

Molnar, S. M., Haiman, Z., Birkinshaw, M., & Mushotzky, R. F. 

2004, Astrophys. J., 601, 22 
Perlmutter, S., et al. 1999, Astrophys. J., 517, 565 
Pierpaoli, E., Borgani, S., Scott, D., & White, M. J. 2003, Mon. 

Not. Roy Astron. Soc, 342, 163 
PLANCK. 2006, the Planck collaboration, "Planck: The scientific 

programme", astro-ph/0604069 
Press, W. H., & Schechter, P. 1974, Astrophys. J., 187, 425 
Rapetti, D., Allen, S. W., & Weller, J. 2005, Mon. Not. Roy. Astron. 

Soc, 360, 555 

Reed, D., Bower, R., Frenk, C, Jenkins, A., & Theuns, T. 2006, 

submitted to Mon. Not. Roy. Astron. Soc, astro-ph/0607150 
Riess, A. G., et al. 1998, Astron. J., 116, 1009 

Robinson, J., & Baker, J. E. 2000, Mon. Not. Roy. Astron. Soc, 
311, 781 

Robinson, J., Gawiser, E., & Silk, J. 1998, astro-ph/9805181 
— . 2000, Astrophys. J., 532, 1 

Ruhl, J. E., et al. 2004, the SPT collaboration, astro-ph/0411122 
Sadeh, S., Rephaeli, Y., & Silk, J. 2006, Mon. Not. Roy. Astron. 
Soc, 368, 1583 

Sasaki, M., Valiviita, J., & Wands, D. 2006, astro-ph/0607627 
Scoccimarro, R., Sefusatti, E., & Zaldarriaga, M. 2004, Phys. Rev., 
D69, 103513 

Seljak, U. 2002, Mon. Not. Roy. Astron. Soc, 337, 769 
Seljak, U., & Zaldarriaga, M. 1996, Astrophys. J., 469, 437 
Spergel, D. N., et al. 2006, submitted to Astrophys. J. Suppl., astro- 
ph/0603449 

Verde, L., Wang, L.-M., Heavens, A., & Kamionkowski, M. 2000, 
Mon. Not. Roy. Astron. Soc, 313, L141 



11 



Wang, S., Khoury, J., Haiman, Z., & May, M. 2004, Phys. Rev., Willick, J. A. 2000, Astrophys. J., 530, 80 

D70, 123008 Zaldarriaga, M. 2004, Phys. Rev., D69, 043508 

Warren, M. S., Abazajian, K., Holz, D. E., & Teodoro, L. 2006, 
Astrophys. J., 646, 881 

White, M. J. 1999, Mon. Not. Roy. Astron. Soc, 310, 511 



